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PREFACE 


This report was prepared for Contract NAS 1-12008, Expansion 
and Extension of the ODIN/RLV Computer Program - Task 2, 
Evaluate and Improve the Existing ODIN Program Library. The 
contract was funded by the National Aeronautics and Space 
Administration, Langley Research Center, Space Systems Divis- 
ion. 

The ODIN procedure is a design analysis technique which allows 
the use of existing computer codes as part of a larger simula- 
tion. Communication of information among computer codes is 
accomplished by means of a data base repository accessible and 
managed by the ODIN executive computer code, DIALOG. 

The objective of the contract task was the improvement of the 
VAMP program and incorporation into the ODIN library. The 
original development on the VAMP computer code was performed 
by LTV Aerospace Corporation under contract to NASA. 
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by 
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and 

C. R. Glatt - Aerophysics Research Corporation 


SUMMARY 

This document describes a method and computer program for 
the calculation of the mass properties , c.g. location, 
enclosed volume, wetted area and planform area of any 
closed structure that has a plane of symmetry. The vehicle 
is described to the computer program by ordered sets of 
X, Y, Z coordinates of points on its surface. The X, Y • Z 
coordinates are converted to quadrilateral elements for 
analysis. The mass properties of each quadrilateral may 
be computed from a thickness and density input for each 
quadrilateral or from a weight per unit area input at each 
point or from a combination of both. The weight per unit 
area can be a composite of the real wall including skin, 
insulation, ribs, stringers, standoffs, brackets, etc. 

The mass properties of the quadrilateral elements are 
accumulated for each section of the vehicle and for the 
total vehicle. The computed mass property totals may 
contain not only the contribution from the distributed mass 
on the vehicle surface wall, but additional masses may be 
added as "black boxes" by specifying each one's c.g. loca- 
tion and mass properties. The added black boxes may be 
inside or outside the surface and do not have to be 
symmetrical with respect to X-Z plane. 

Program VAMP (Volume, Area and Mass Properties) can produce 
picture-like images of the vehicle or individual sections from 
the quadrilateral element data. This facilitates checking the 
input and visualizing new/or modified configurations. Ortho- 
graphic, perspective and stereo views may be obtained. 

The VAMP program has been integrated into the ODIN (Optimal 
Design Integration) Library of programs so that any selected 
mass properties data can be stored in a common data base for 
use by other programs. Alternately, VAMP may accept data 
from the common data base as input to the program. 



INTRODUCTION 


The estimation of the mass properties of a vehicle is one 
of the most important considerations in the design process 
and yet one of the most inexact engineering endeavors. 

While the calculation of aerodynamic, propulsion and 
mission performance are based on widely recognized mathe- 
matical prediction techniques, the estimation of weight 
must be based largely on historical data. The art of 
weight estimation has evolved through the years by the 
diligent collection and correlation of component weights 
of previously built vehicles. New design weights are 
predicted on the basis of the component weights of past 
designs. Little information is usually available on the 
other properties such as volume, area, center of gravity 
and inertia of the components. Even for "weights only" 
calculation, historical methods lack validity in predict- 
ing advanced state-of-the-art design data. Further, the 
influence of small design perturbations eludes such gross 
prediction methods. 

During the later design phases when detail structured 
analyses are available for some components, the weights 
engineer must still revert to empirical analysis methods 
for subsystem and secondary structure weights such as 
wing leading edges, canopies, fairings, etc. 

Existing programs such as references 1 through 6 lean 
heavily on the historical approach for weight estimating. 

To obtain more reliable information, the designer must 
turn to structural analysis and subsystem design programs 
such as references 7 through 13. The detailed geometry 
descriptions required for these classes of analyses are 
usually not known until the later design phases. At this 
later phase the design is usually frozen and little 
geometric perturbation is allowed. Further, small per- 
turbation analysis using these programs is very expensive 
in terms of manpower and computer resources. 

Although the precise weight and other properties of the 
vehicle components are illusive, or at best a very complex 
calculation, mathematical equations for combining the com- 
ponents .mass properties into total vehicle mass properties 
are quite concise. Mass property evaluations for the 
total vehicle as well as perturbations due to added or 
deleted masses are easily determined once the mass properties 
of all the components are known (or assigned) . 
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The program VAMP (for Volume Area and Mass Properties) is 
designed to bridge the gap between the historical weight 
approach and the detailed structural analysis approach. 
VAMP computes the volume, area and mass properties based 
on a surface model and a surface unit distribution of 
mass. Any number of surface models or sections may be 
described. The program maintains an accounting of the 
properties of each component and can accumulate the 
properties of all vehicle sections. The accumulated 
section properties can be combined with the properties 
of arbitrary specified masses (black boxes) . 

The original development on the VAMP computer code was 
performed by Mr. Patrick Norton under contract to NASA. 
Robert T. Wingate served as contract monitor. Mr. Phillip 
J. Klich, also of NASA, made significant contributions to 
the development and final checkout of the VAMP computer 
program. VAMP was later modified for use in the ODIN 
system and has since been used extensively in design 
analysis studies at NASA/Langley Research Center. 


The significant contribution of the VAMP program develop- 
ment to the general problem of computer aided design is 
twofold. First, the program offers a unified approach to 
the accounting of the vehicle component mass properties. 
The fact that a component weight is not well known does 
not exclude it from being analyzed in a mass properties 
evaluation. Furthermore, detailed mass properties of 
some vehicle components may be determined elsewhere and 
still be included in the VAMP analysis. Second, the 
distributed mass surface model incorporated in VAMP 
allows a very convenient means of producing vehicle 
geometric perturbations early in the design process with 
some degree of confidence in the relative affect of the 
perturbation. As more detailed information on the design 
becomes available, the component input to VAMP may be 
gradually replaced without distruction of the overall 
continuity of the weights analysis. 
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VOLUME, AREA AND MASS PROPERTIES ANALYSIS 

The mathematical modeling of the VAMP program permits the 
calculation of the mass properties, c.g. locations, 
enclosed volume, wetted area and planform area of any 
closed structure that has a plane of symmetry, e.g. fuse- 
lage, stiffened fuel tank, etc. The vehicle is described 
in a symmetrical manner and consists of sets of X, Y, Z 
coordinates on its surface. The surface coordinates are 
converted to quadrilateral elements but since the four 
corners of the quadilaterals are not necessarily coplanar, 
each quadrilateral is analyzed as two triangles. The mass 
properties are computed from a thickness and density input 
for each quadrilateral, from a weight per unit area input 
at each point, or from a combination of both. The weight 
per unit area can be a composite of the real wall includ- 
ing skin, insulation, ribs, stringers, standoffs, brackets, 
etc . 

The elemental mass properties are accumulated for each 
section and for the total vehicle. The computed mass 
property totals containing the contribution from the 
distributed mass in the vehicle surface wall may be 
combined with additional masses specifying each one's 
c.g. location and mass properties about its c.g. The added 
black boxes may be inside or outside the surface and do 
not have to symmetrical. 

The properties of the vehicle computed by the program are: 

1. Weight 

2. Coordinates of the c.g. 

3. Three mass moments of inertia plus I v „ about 

both the reference axes and axes through the 
c.g. parallel to the reference axes. 

4. Surface area. 

5. Enclosed volume. 

6. Coordinates of the c.v. 

7. Projected area. 

These properties are output for each segment along with 
cumulative totals and complete vehicle values. 
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Surface Model 


The numerical model of the surface shape to be analyzed is 
described by. quadrilateral elements. Each quadrilateral 
element consists of a grouping of four surface points. An 
organized set of quadrilateral elements form a segment of 
the configuration. A number of segments may be used to 
provide a complete description of the configuration. 

The configuration is usually (but not necessarily) positioned 
with its nose at the coordinate system origin and with the 
length of the body stretching in a positive X-direction. 

The coordinate system used for this program is the right 
hand Cartesian system as illustrated below: 


Z 



The configuration is assumed symmetrical about the X-Z 
plane, and the analysis accounts for the plane of symmetry. 

The configuration may consist of an arbitrary number of 
wing and body sections. The body sections are described 
in the positive half plane perpendicular to the plane of 
symmetry. The wings are described as airfoil sections 
parallel to the plane of symmetry. Canards and horizontal 
stabilizers are described similar to wings. Fins cannot 
be described as airfoil sections but they may be described 
as body sections. 
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Points along the X-axis are referred to as station coordinates 
A Y-Z plane is a station plane . The Y-Z coordinates in the 
station plane form a station contour. 


-tz 



A portion of the mathematical surface model is presented 
in the illustration above. The X,Y,Z coordinates of the 
points where the solid grid lines cross are described to 
the program. Each point (I, J) is associated with a station 
number J and a contour point I. Thus, the circled point 
is 8,5 and it has coordinates X,.; Yg , Z g , 5 . The model 
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groups the surface points into quadrilateral elements which 
also have unique I,J identification. The crosshatched one 
is 8,5 corresponding to its lowest numbered corner. The 
dashed diagonals illustrate the division of each quadri- 
lateral into two triangles for analysis. 

Airfoil sections are described in a similar manner as the 
body sections above. In the case of airfoil section, the 
parts are described in X-Z planes parallel to the plane of 
symmetry. 

The objective in the analysis is to calculate the properties 
of each quadrilateral (two triangles) elements and trans- 
form the results back to the reference coordinate system, 
then accumulate the properties from all elements of the 
surface. 

Mass Properties of a Triangular Surface Element . - The 
triangular surface element illustrated below with its base 
a station J and apex at J+l is assumed to be of uniform 
thickness and homogeneous. 


y 



As viewed from the outside of the vehicle, the local x,y,z 
coordinate system has its origin at (I,J) . The x-y plane 
is coincidental with the plane of the triangle. A is the 
height of the triangle. B is the base of the triangle. 
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C is the offset of the apex from the y-axis. H is the 
thickness which is assumed to be small compared to A, B 

and C. 

The area, S of this elemental triangle is: 


S = AB/2 (1) 

The first moments of inertia for the same triangle are: 
S = PHS (B+C) /Z (2) 

X 

S y = PHSA/2 (3) 

S z = 0 (4) 


The second moments of inertia for the elemental triangle 

are; 


*xx - PHA 2 S/6 

(5) 

Iyy = pH(B 2 +BC+C 2 )S/6 

(6) 

I =1+1 

z z xx yy 

(7) 

Where p is the density of the thin triangle 
products of inertia are: 

. Finally, 

I = pH (B+2C) AS/12 

xy 

(8) 

I =0 

yz 

(9) 

I =0 

xz 

(10) 

The coordinates of the centroid (center of 
triangular element are: 

gravity) of 


x = (B+C)/3 


( 11 ) 



The transformation of the inertias and products of inertia 


to the centroid 
I 

XX 

yield: 

SA 2 /18 

(13) 

T 

X (B 2 -BC+C 2 )/18 

(14) 

ii 

>1 >i 

>i , X 
|H 

SA(ZC-B)/36 

(15) 

I 

I +1 

(16) 

ZZ 

xx yy 


Transformation Matrix for Triangle Elements , - The mass 
properties above are given for a simple triangular element 
in the local x,y,z coordinate system. The elemental mass 
properties must be transformed to the reference coordinate 
system X,Y,Z for cummulation with other elements. The 
transformation from the primary X,Y,Z coordinate system to 
the local system is given by: 


X 


"x - xo" 

y 

i — i 

a 

uEi 

ii 

Y - YO 

_ z „ 


_z - zo. 


Where XO,YO,ZO are the coordinates of the prime system 
origin with respect to same arbitrary reference point and 
DC is a matrix of direction cosines. The elements of DC 
are : 



XPX 

XPY 

XPZ" 

II 

YPX 

YPY 

YPZ 


_ZPX 

ZPY 

ZPZ. 


For a unit vector in the x direction, XPX is the projection 
of x on X , XPY is the projector of x on Y , etc. The 
vertices of the triangular element in local x,y,z 
coordinates may be considered as vectices B and D, drawn 
from the origin (0,0,0). The components of B and D are 
developed as follows. That is,BX is the X component of 
B and DY is the Y component of D, etc. Therefore: 
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If: XO = X(J) (19) 

YO = Y (I , J) (20) 

ZO = Z (I, J) (21) 

Then: BX = 0.0 (B lies in YZ Plane) (22) 

BY = Y (1+1 , J) - YO (23) 

BZ = Z(I+1,J) - ZO (24) 

DX = X(J+1) - XO (25) 

DY = Y (1+1 , J+l) - YO (26) 

DZ = Z (I-1,J-1) - ZO (27) 

The cross product of B with D is a vector S2 in the z 
direction. Its magnitude is twice the surface area of 
the triangle. Let SX,SY and SZ be the components of S2. 
Since BX = 0, the cross product yields: 

SX = BY*DZ - DY*BZ (28) 

SY = BZ*DX (29) 

SZ = -DX*BY (30) 

Let B, D, S2 represent the magnitude of their respective 
vectors. Then IT/B = unit vector in x direction and S2/S2 
= unit vector in z direction. 

It follows: 

S2/S2 cross B/B = unit vector in 

y direction (31) 

The components of each unit vector are a row in the matrix 
DC. Therefore: 

XPX = BX/B =0.0 (32) 

XPY = BY/B (33) 

XPZ = BZ/B (34) 

ZPX = SX/S2 (35) 

ZPY = SY/S2 (36) 

ZPZ = SZ/S2 (37) 

YPX = ZPX*XPZ - XPY* ZPZ = DX/A (38) 
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YPY = XPZ*ZPX 
YPZ = ZPX*XPY 


(39) 

(40) 


The preceding steps considered a triangle with its base on 
J and apex on J+l. The other half of each quadrilateral 
has its base on J+l and vertex on J. The vertices are 
identified as follows: 


(I, J+l) 



The analysis for the inverted triangle is similar to the 
upright triangle. The only difference is in equations ( 13 ) 
through (21) . For the inverted triangle. 


XO = X(J+1) (41) 

YO = Y (1+1 , J+l) (42) 

ZO = Z (1+1, J+l) (43) 

BX = 0.0 (44) 

BY = Y (I, J+l) - YO (45) 

BZ = Z (I, J+l) - ZO (46) 

DX = X(J) - XO (47) 

DY = Y (I , J) - YO (48) 

DZ = Z (I, J) - ZO (49) 


The above analysis is similar for airfoil type sections 
described in the x-z plane. 
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Mass Properties of the Surface. - The c.g. and mass pro- 
perties expressions for a triangle are given in the 
previous section. The triangular element parameters A 
and B have been developed above and are always positive. 

2 2 2 

Now C =■ D - A but C can be negative. However, C and 
its sign are given directly by the dot product of D with 
the unit vector in the x direction. 

C = DY*XPY+DZ*XPZ (50) 

The resulting c.g. coordinates are in the x, y, z system 
and the inertias are about axes through the c.g., parallel 
to the prime system. They are transformed to the system 
reference axes as follows: 


xbr" 


xo" 


X 




-1 



YBR 

— 

YO 

+ [DC] * 

y 

_ZBR 


_Z° 


_ z 


This is readily expanded since the inverse of DC is 
identical to its transpose. The products and moments of 
inertia parallel to the prime axes are assembled in PMIP: 



X 

X 

H 

1 

"" I XY 

-I xz 

PMIP = 

X XY 

X 

X 

H 

1 

_I YX 


- J xz 

_I YZ 

_I zz- 


The minus signs on the products of inertia are a consequence 
of the coordinate system definition. This places the matrix 
PMIP In tensor form. It can be shown that the inertias are 
rotated parallel to the system axes by: 

[PMI ] = [DC] T * [PMIP]* [DC] (53) 

PMI contains the inertias about axes parallel to the system 
reference axes, but centered on the triangles c.g. Trans- 
lation to the system origin makes use of the parallel axis 
theorem along with XBR, YBR, ZBR computed above, and is 
included with the summation procedure described in the next 
paragraph. The products of inertia contained in PMI are 
the negative of their engineering definition, like they 
are in PMIP. 
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In the VAMP analysis, all values in PMI are twice the 
true value since S2 = twice the area of the triangle, S2 
has been retained through all the computations. But each 
triangle has a symmetrical counterpart on the other side 
of Y = 0 whose contribution to the total is either equal 
or equal- and-opposite. In the latter case, the resultant 
is zero; i.e., I XY = I = 0.0. 

Area Calculations . - The surface area and projected area 
result from the summation of the elemental areas from the 
above analyses. The surface area is: 

A = 1/2 2 S2 i (54) 

S2 is twice the elemental area and the projected area (on 
the y-z plane) is: 

A = 1/2 2 SY. (55) 

P i 1 

Volume Calculations . - The elemental volume of a surface 
element is the projection of the elemental area on the 
x-y plane multiplied by the average distance to the c.g. 
of the element. The total volume is: 

V = 2 x. • ZPX. • A. (56) 

i 

The sign of the direction cosine accounts for whether 
the elemental volume is added or subtracted. 

The center of volume is simply: 

CV = ^ Z X 2 . ZPX. A. 

V l ll 

The above calculations are approximate because of the 
assumption that the c.g. of the element is an adequate 
approximation for the entire elemental surface. The 
assumption is reasonable for small elements but may 
not be valid for a gross panelling of the configuration. 
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Black Box Contributions to the Mass Properties 


The surface model described above is suitable for describ- 
ing shell structures such as tanks, skins, etc. However, 
the addition of mass not conforming the above model may be 
necessary in the complete description of the vehicle mass 
properties. In the VAMP analysis, additional mass sources 
may be added by specifying each ones center of gravity 
(c.g.) location and mass properties. These mass sources 
called "black boxes" may lie inside or outside the sur- 
face and do not have to be symmetrical with respect to 
the x-z plane. The VAMP program combines the detailed 
shape inputs and the black box data to produce the over- 
all mass properties of the vehicle. 

Black box masses are described to the program in the 
following manner: 

1. Location of its c.g. in system reference 
coordinates. 

2 . Total weight. 

3 . Moments and products of inertia about local 
axes through its c.g. 

4 . Orientation of local axes relative to the 
system reference axes or relative to some 
intermediate axes. 

The minimum required for a black box input is the specifica- 
tion of items 1 and 2. If items 3 and 4 are omitted, 
the moments and products of inertia are assumed to be 
zero and the orientation of the local axis is assumed 
to be coincidental with the system reference axes. The 
contribution of each black box is added to the cumulative 
totals produced by the analysis. Since the size and shape 
is not being specified, there is no addition to surface 
area or volume. Since the black boxes do not have to be 
symmetrical with respect to the x-z plane, I^ Y ^ I YZ ^ 

0.0, and there can be a non-zero first moment with respect 

to Y . 

The moments and products of inertia of each black box are 
input directly to PMIP (see equation 52) . The program 
takes care of filling the upper half and adding the minus 
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signs . Rotation of the inertias uses equation 53 . DC 
is obtained from input data as follows. Consider the 
local x,y plane: 



Point £ is the origin. Points 1 and 2 are arbitrary 
points on the +x and +y axes respectively. Coordinates 
of these three points, measured in the system to which 
the inertias are being rotated, are inputs. It is a 
simple matter to construct unit vectors in the x and y 
directions and their components yield the first two rows 
in DC , equation 18. The third row is the cross product 
of these two unit vectors. 

In order to construct DC, it is not necessary that point 
0 be the c.g. . The coordinates x , y must be parallel to 
the axes to which the PMIP inputs are referenced, but 
rotation is performed independently of translation. After 
the rotation of equation 53, the inertias in PMI are still 
centered on the c.g. of the black box. Only the rotation 
takes place and the contents of PMI are moved to PMIP . 
Thus, it is possible to make several rotations, if needed, 
to reach alignment with the system reference axes. The 
input to the program can be developed from any convenient 
points on the black box and intermediate structure whose 
coordinates are known. The input does not always have 
to be in the same units as the system X , Y ,Z data. 
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Configuration Plotting 


The plotting capability in the VAMP program is based 
upon reference 16. The following types of plotting 
capability are contained in the program: 

1. Three-views. 

2. Orthographic, from an arbitrary viewing angle. 

3. Perspective, from ah arbitrary viewing angle. 

4. Stereoscopic, from an arbitrary viewing angle. 

The program interfaces through the CDC 6600 to the follow- 
ing types of equipment: 

1. On-line cathode ray tube. 

2. CALCOMP plotter. 

3. Gerber plotter. 

4. Stereoscope. 

The numerical model of the aircraft configuration may 
include any combination of components available in the 
basic VAMP program. The wing is made up of airfoil 
sections , the body is defined by either circular or 
arbitrary sections. The fins must be defined similar 
to the fuselage, but stabilizers and canards can be 
defined similar to the wings. The vehicle geometric 
specification is converted into a set of quadrilateral 
panel elements in a manner similar to that described 
above. However, quadrilateral elements are analyzed 
directly rather than as two triangles. 

Orthographic projections are created by rotating each 
point on the body surface to the desired viewing angle 
and then transforming the points into a coordinate system 
in the plane of the paper. The rotations of the body and 
its coordinate system to give a desired viewing angle are 
specified by angles of roll, pitch and yaw (cj>, 6, \jj) using 
the convention below: 
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The code computes the "average" unit normal vector to each 
panel. The resulting set of vectors may be used to pro- 
vide the capability of deleting most elements on the sur- 
face of the configuration which would not be seen by a 
viewer. By this device a user may remove many confusing 
panel elements. No provision is made for deleting com- 
ponents hidden by other components or for deleting portions 
of an element at the present time. 

When three-views are requested, the plan, front and side 
views are provided in a compact and pleasing to the eye 
arrangement. An option is provided for the orthographic 
projections of these three-views to be spaced one above 
the other. A typical three-view obtained in this manner 
has been presented in figure 1. 

The perspective views represent the projection of a given 
three-dimensional array. The two-dimensional view is 
constructed relative to a viewing point and a focal point 
specified by coordinate points in the data coordinate 
system. Data are scaled to the viewer page size auto- 
matically by the specification of the viewing field 
diameter and the viewing field distance. The coordinates 
of the viewing point determine the position from which 
the data array will be viewed and the coordinate values 
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TYPICAL THREE-VIEW FROM VAMP, 



of the focal point control the direction and focus. The 
size of the projection on the viewing plane reflects the 
distance between the viewing point and the focal point. 

Data which are within the cone of the viewing plane but 
not in the immediate range of the focal point may be 
distorted. Perspective may be eliminated by specifying 
a large viewing field distance. A typical detailed ortho- 
graphic projection of a shuttle type vehicle is presented 
in figure 2. 

The above explanation of the perspective plots also applies 
to the stereo views. The use of the stereo option causes 
the program to be executed twice in setting up two plots 
for the left and right frames. These frames are suitable 
for viewing in a stereoscope. 
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FIGURE 2 TYPICAL ORTHOGRAPHIC PROJECTION WITH HIDDEN LINES INCLUDED, 


PROGRAM USAGE 


The VAMP input data is organized by cases as illustrated in 
figure 3. Each case consists of a title card, configuration 
shape data, black box input data and plot data. Either or 
both configuration shape data and black box data may be 
specified. Plot data is optional and its use depends on the 
values of certain input parameters. If both configuration 
and black box data are specified, the black box data follows 
the shape data. Plot data is always at the end of the input 
list. 


Control Cards 

The program case data is preceded by the control cards 
required to retrieve the VAMP program from permanent stor- 
age and execute the program. Figure 4 illustrates the con- 
trol cards to use VAMP at Langley Research Center (LRC) . 
Figure 4A illustrates execution from a stored machine 
language program. Figure 4B illustrates a compile, load 
and execute run from stored source code assuming program 
modifications are required. Figure 4C illustrates the use 
of VAMP within the ODIN system. In the ODIN system VAMP 
is stored in absolute binary form. 

The dashes on the JOB, USER and REQUEST cards indicate other 
information which the user must supply in accordance with 
LRC computer complex accounting procedures. Further, the 
wedge number (program storage location) is subject to change. 
The latest wedge number is available from Langley (LRC) 

The user must supply the estimated run time in CPU seconds, 
field length for the job and 0/S call count on the JOB card. 


Typical values for a 

single 

case are tabulated 

below: 


CPU 

FIELD LENGTH 

O/S CALLS 

EXECUTE ABSOLUTE 
BINARY PROGRAM 

20 

42000 

500 

COMPILE, LOAD 
AND EXECUTE 

40 

51000 

800 

EXECUTE WITHIN 
ODIN SYSTEM 

60 

56000 

1000 


The tabulated values for CPU and 0/S calls estimated in the 
above tables will vary with the complexity of the configura- 
tion data, the number of cases and the number of plots. A 
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FIGURE 3 ILLUSTRATION OF VAMP DECK SETUP. 



JOB ,1,20*41000,500. 

USER. - 

FETCH ,A36 8 3 , SPR , BINARY, ,OVAMP 2 . 

0 VAMP 2. 

REQUEST, TAPE9 8, HI. CALTP,RIL, 

REWIND f CALTPE,TAPE99. REQUIRED ONLY IF 

COPYBF , CALTPE, TAPE99 . PLOTS ARE REQUESTED 

UNLOAD, TAPE 9 9 . 

7-8-9 

(VAMP CASE DATA) 

S-7-8-9 

FIGURE 4A EXECUTION OF STORED PROGRAM 


JOB, 1,40, 51000, 800. 

USER. 

FETCH , A3 6 8 3 , SPR , SOURCE . 

RUN , S , , ,SCFILE. 

LGO. 

REQUEST, TAPE9 8, HI. CALTP ,RIL, 

REWIND, CALTPE, TAPE9 9. REQUIRED ONLY IF 

COPYBF, CALTPE, TAPE9 9. PLOTS ARE REQUESTED 

UNLOAD, TAPE 9 9. 

7-8-9 

(MODS TO SOURCE PROGRAM, IF ANY) 

7-8-9 

(VAMP CASE DATA) 

6-7-8-9 

FIGURE 4B COMPILE, LOAD AND EXECUTE 


'EXECUTE VAMP' 
(VAMP CASE DATA) 

7-8-9 

'EXECUTE PLOTS V' 
7-8-9 


FIGURE 4C EXECUTE VAMP WITHIN ODIN SIMULATION 


FIGURE 4 ILLUSTRATIONS OF CONTROL CARDS REQUIRED FOR VAMP. 


REQUIRED ONLY IF 
PLOTS ARE REQUESTED 
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good rule to follow is to allow ample CPU and 0/S calls on 
the first run, then use the results to estimate these para- 
meters of future runs. 


General Input Procedure 

A typical VAMP input deck will consist of four types of input 

1. Title card in free-field format. 

2. Surface model data ($SHAPE) in NAMELIST format. 

3. Black box input data ($BLKBOX) in NAMELIST format. 

4. Plotting information in fixed field format. 

These input sets may be combined in various ways depending 
upon the input parameters in $ SHAPE and $BLKBOX input sets. 
All input information, including NAMELIST, must be placed 
in the first 72 columns of the input cards. Columns 73 
through 80 are reserved for comments and identifiers. The 
purpose of this section is to summarize the input data 
requirements and describe the input data flow logic. 

Title Card - The title card is a single card with any desired 
information in columns 1 through 72. The contents of this 
card are printed at the top of each page of output and on all 
plots. 

$SHAPE Input Set - The $SHAPE input data description in 
figure 5 contains the corner-point geometry of the component 
external surfaces, the elemental surface thicknesses, the 
surface weight distribution, etc. The configuration is 
described by ordered sets of points on the surface for each 
component. Each component is subdivided into sections for 
input purposes. Any number of components may be described 
to the program. All of this surface model data can be 
supplied in the $SHAPE input list or some may be supplied 
by Harris type input (reference 16) or digitizer input as 
described later. A minimum of two $SHAPE input sets is 
required for each segment. If any of the auxiliary data 
options, KXYZ, KRHOH or KSWT are invoked they must be set 
in the first $SHAPE set and the appropriate supplemental 
formatted data cards must be placed between the first and 
second $SHAPE set. The second $SHAPE set normally terminates 
the input for one component of the configuration. However, 
additional $SHAPE sets may be required if the user decides 
to rearrange the section data (KEEP option) . 
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NAME 

DEFAULT 

VALUE 


DESCRIPTION 

CNVMS 

32.175 

(A) 

Factor to convert weight to mass. 

H (40 , 25) 

1000*0. 

(D) 

Surface element thickness. 

KEEP (30) 

30*0 

(C) 

Sequence of contain points to be used- 
in constructing a component. 

KPLT 

0 

(A) 

Compute/plot central indicator. 

KPLT = 0; computes mass properties 
and plots. 

KPLT =-l; computes mass properties only. 
KPLT =+l; plots only. 

KRHOH 

0 

(A) 

Density input option if KRHOH = 1, 

RHO and H are input in 10F7.0 format, 
after $SHAPE input list. 

KSWT 

0 

(A) 

Surface weight input option if KSWT = 1, 
elemental surface weights are input 
in 10F7.0 format after $SHAPE input list. 

KWNG 

0 


Wing input option (10F7.0 format) 

KWNG = 0; no wings. 

KWNG = lj define upper and lower surface; 
KWNG =2? define upper surface and 
mean camber. 

KWNG = 3 . define both surfaces. 

KXYZ 

0 (A 

rE,) 

Corner point geometry input option. 

KWYZ =0; data expected in $SHAPE list. 
KXYZ = 1; data expected in 10F7.0 format. 
KWYZ = 2 ; data for circular cross section. 
KXYZ = 3; digitizer input. 

KYZRO 

0 

(B) 

Contour closure option. 

KYZRO = 0; program closes contour to 
x-y plane. 

KYZRO = 1; program sets last y,z value 
to the first y,z values. 
KYZRO = 2; leaves the input unchanged. 

LABEL (6) 

0 

(B) 

Label of 6 words for each segment. 

NOSE 

0 

(B) 

Indicator identifying a singular station 
at the nose (constant y,z - values). 

NPTS 

0 

(A) 

Number of points on each station contour 
(2<NPTS<45) . 


FIGURE 5A $SHAPE INPUT LIST, 
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NAME 

DEFAULT 

VALUE 


DESCRIPTION 

NS EG 

0 

(A) 

Number of segments or components for 
the configuration. 

NSTA 

0 

(A) 

Number of stations in the current seg- 
ment or component (2>NSTA>25) . 

RFLNTH 

100. 

(A) 

Reference length for the vehicle or 
configuration. 

RHO (40,25) 

1000*0 

(D) 

Elemental density. 

SCLX 

1.0 

(A) 

X-scale factor for geometric input. 

SCLY 

1.0 

(A) 

Y-scale factor for geometric input. 

SCLZ 

1.0 

(A) 

Z-scale factor for geometric input. 

SURFWT (40,25) 

1000*0 

(D) 

Elemental units weights for the 
surface components. 

X (40) 

40*0 

(D) 

X-position of body station planes. 

Y (40,25) 

1000*0 

(D) 

Y-position of body station planes. 

Z (40,25) 

1000*0 

(D) 

Z-position of body stations. 

XOF 

0 . 

(A) 

X-offset from system reference. 

YOF 

0 . 

(A) 

Y-offset from system reference. 

ZOF 

0 . 

(A) 

Z-offset of system reference. 

XWNG (40,25) 

1000*0 

(D) 

X-position of wing station planes. 

YWNG (40) 

40*0 

(D) 

Y-position of wing sections. 

ZWNG (40,25) 

1000*0 

(D) 

Z-position of wing section. 

A. Reset 

before each 

new 

case. 


B. Reset before each segment. 

C. There is a second $ SHAPE set required for every segment. 
KEEP may be input in this second set. If KEEP is non-zero, 
the program returns and requests another "second set." 

KEEP is reset to zero prior to reading each second set. 

D. Values can be indefinite depending on the loader option 
set. Indefinite means the initial value will be the value 
in the computer after the program is loaded or after a 
change in overlays during execution. 

E. The Y,Z arrays generated when KXYZ = 2, can be retained 
for the next segment by changing to KXYZ = 0 in the next 
segment. 


FIGURE 5B JSHAPE INPUT LIST (CONTINUED), 
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The arrangement of the $SHAPE data is illustrated below: 

$SHAPE -$ 

Optional formatted data 
$SHAPE- - - -$ 

Optional rearrangement of sections 
$SHAPE -$ 

$BLKBOX Input Set . - The $BLKBOX input contains the informa- 
tion necessary to describe auxiliary masses which cannot be 
described by geometric input data. Any number of black boxes 
may be described to the program and they may be symmetrical 
or nonsymmetrical with respect to the configuration center 
plane. Figure 6 briefly describes the input parameters in 
the $BLKBOX input. 

Plotting Information . - Four CALCOMP plot options are avail- 
able in VAMP. 

1. Orthographies projection. 

2. Three views. 

3. Prospectives . 

4. Stereo frames. 

Each of the above options may be invoked by a single formatted 
input to be described later. The components plotted are those 
specified in the $SHAPE input sets above. 


Input Flow Logic 


Figure 7 shows the input flow logic for the program. After 
initializing the data, VAMP reads a title card, then proceeds 
to read in the surface model (figure 7B) . The first $SHAPE 
data set determines the flow path within the VAMP program for 
the component. If the number of components (NSEG) is less 
than one, program flow returns (figure 7A) assuming there is 
no external geometric configuration data. If one or more 
segments are to be input, the auxiliary data option parameter 
(KXYZ) is tested. Three auxiliary data options can be invoked 
depending upon the value of KXYZ. 

1. Read Harris-type input (KXYZ = 1). 

2. Compute circular cross section (KXYZ = 2) . 
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NAME 

DEFAULT 

VALUE 

DESCRIPTION 

BXPLN (3 , 3) 

9*0. 

(C) 

Coordinates of three points. 

CNVMS 

32.174 (A) 

Factor to convert weight to mass. 

IXX 

0 . 

(B) 

Black box inertia about load x-axis. 

IYY 

0 . 

(B) 

Black box inertia about load y-axis. 

IZZ 

0 . 

(B) 

Black box inertia about load z-axis. 

LABEL ( 6 ) 

0 . 

(B) 

Alphanumerical label up to 6 words. 

NBX 

0 . 

(C) 

Controls transformation of black 
box to reference axis. 

NBX = 0; no (more) black boxes. 

NBX = 1; transform inertias (PMIP) 
to reference system and 
accumulate . 

NBX = 2; If black box has symmetrical 
counterpart. 

NBX =-l; Rotate inertia but not 
translated. 

PMIP (3,3) 

9*0. 

(B) 

Black box moments and products of 
inertia . 

WTM 

0 . 

(B) 

Black box weight or mass 

xo 

0 . 

(a,d) 

X-position of black box c.g. 

YO 

0 . 

(A, D) 

Y-position of black box c.g. 

zo 

0 . 

(A, D) 

Z-position of black box c.g. 

XOF 

0 . 

(A) 

X-offset from the system reference 
axis . 

YOF 

0 . 

(A) 

Y-offset from the system reference 
axis . 

ZOF 

0 . 

(A) 

Z-offset from the system reference 
axis . 


A. Reset before each new case. Thus, the value is carried 
forward from the vehicle wall computations. 

B. Reset before each new black box. 

C. Reset before each $BLKBOX set is read. NBX must appear 
in every $BLKBOX set except the last one. 

D. Only the second and third columns of BXPLN are reset 
to zero. The first column contains XO,YO,ZO which 
retain their input values. They are not changed by 
the offsets (XOF, etc.) 

FIGURE 6 $BLKB0X INPUT LIST. 
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FIGURE 7 A INPUT FLOW LOGIC. 
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FIGURE 7B INPUT FLOW LOGIC. 
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3. Digitalized data cards (KXYZ = 3) . 

The first and third options require additional input which 
automatically overrides the $SHAPE input geometry. Density 
and thickness may also be read in the auxiliary data format 
depending on the value of KRHOH. Surface weight data may 
be read by the auxiliary data format by use of the input 
parameter KSWT. The stacking of the auxiliary data described 
above can be determined from the input flow logic (figure 7B) . 
In general, the sequence is: 

1. Auxiliary geometry data. (KXYZ>0) 

2. Auxiliary density/thickness data. (KRHOH>0) 

3. Auxiliary surface weight data. (KSWT>0) 

If any of the auxiliary data options are not invoked, the 
data subset is simply omitted. 

After all auxiliary data is read into the program, a second 
$SHAPE input set is read by the program. The second set pro- 
vides the user the opportunity to rearrange the sections of 
the current vehicle component. This is done through the use 
of the KEEP option. If the KEEP option is invoked, the 
sections are rearranged according to the sequence provided 
in the KEEP array, then another $SHAPE set is read by the 
program. After all KEEP options are satisfied, the program 
proceeds to add offsets and scale the corner-point geometry. 
Then the data is modified by the closure option, KYZRO. If 
the KPLT option is invoked, the data for the current section 
is written on an auxiliary file for later use by the plot 
subprogram. 

After completion of the above sequence, the program proceeds 
to check for $BLKBOX data. Any number of black boxes may be 
added through multiple $BLKBOX data sets. After all black 
box data is read and accumulated with the rest of the vehicle 
mass property data, the program proceeds to the plot option. 

If plot cards are present, the program plots the information 
contained on the auxiliary plot file. After all the plot 
request cards are satisfied, the program proceeds back to 
the initialization of the data and a new case. If no case 
data is present, the program stops. 
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Configuration Shape Data 


The surface of the configuration being analyzed is described 
to the program by ordered sets of corner-point geometry. 
Components can be sectioned for input purposes by one of the 
methods illustrated in figure 8. Body stations (Y-Z plane) 
or wing stations (X-Z plane) may be described depending upon 
the value of the KWNG parameter. In either case, the pro- 
gram assumes symmetry about the X-Y plane so that only one- 
half of the configuration geometry need to be input. For 
the body stations the right half (facing aft) is input. The 
right wing airfoil sections are described to the program. 

The mass properties of the surface are supplied to the pro- 
gram in elemental form. The weights may be input directly 
as elemental surface weight contributions or they may be 
input as density and thickness specifications of each surface 
panel. Elemental mass properties may be specified in NAMELIST 
format or by formatted input cards similar to geometric data. 
Fixed offsets and geometric scale factors may also be supplied 
for each component. 

Input Specifications for the $SHAPE Input Set . - The normal 
mode of input is through the $SHAPE input list which includes 
the input parameters for many of the program options. All 
geometric and mass properties data may be specified through 
this mode of input. No auxiliary data read is essential to 
the operation of the program. 

The use of NAMELIST input was chosen for the following reasons: 

1. It is a simple name oriented input easily understood 
by most computer users. 

2. The format is standard and does not require relearn- 
ing from program to program. 

3. It is easily modified by the engineer or programmer 
when adding input variables to the program. 

When a NAMELIST read is encountered in the program, the 
entire input file is scanned up to an end-of-file or a 
record with a dollar in column 2 followed immediately by 
the NAMELIST name requested by the programs. Suceeding data 
items are read until a second dollar is encountered signify- 
ing the end of the NAMELIST. Any data on the input file 
before the requested NAMELIST is found will be ignored. All 
data between the opening and closing dollar is interpreted 
by the NAMELIST input routine. The data item within the 
NAMELIST statement may be in any of three forms: 
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V = c. 


A — D^,...,Dj, 

A (n) = D. , . . . ,D , 

V is a variable name; C is a constant; A is an array name 
and n is an integer constant subscript, are simple 

constants or repeated constants of the form k*C, where k is 
the repetition factor. Data items and constants must be 
separated by commas. The number of constants, including 
repetitions given for an unscripted array must equal the 
number of elements in that array. For a subscripted array 
name, the number of constants need not be equal but may not 
exceed the number of array elements needed to fill the array. 
More than one card may be used for input data and arrays may 
be split between cards. All except the last record must end 
with a constant followed by a comma and no sequence numbers 
may appear. The first column of each record is ignored. 

The set of data items may consist of any subset of the var- 
iable names associated with the NAMELIST name and the name 
need not be any particular order. More details on the use 
of NAMELIST are available in any FORTRAN users guide, but 
the above description should be sufficient for the operation 
of the VAMP program. 

The following paragraphs provide detailed descriptions of 
each variable name in the $SHAPE input list in alphabetical 
order. The descriptive material represents an expansion of 
the abbreviated description given in figure 5. 

CNVMS : Factor to convert weight units into mass units. This 

input parameter affects only the moment and product of inertia 
calculations. The total weight remains in input units. The 
weight elements entered into an inertia calculation are 
divided by CNVMS for conversion to mass. The initialization 
of CNVMS is preset to 32. 1741b/slug. 

H ( 40,25) : Thickness of the quadrilateral surface elements. 

T(T7j) Is the thickness of an element defined by four points 
starting at the geometric point (I,J) and extending one 
point beyond in each direction. It will require change only 
if units of mass other than English units are used. Any 
quadrilateral that has H<_0.0 is deleted from the surface area 
and mass properties calculations, but the volume contribution 
is accumulated for all elements. However, there must be at 
least one non- zero thickness between every station if .mass 
properties are to be computed. 
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KEEP (point) ; An array of integers defining a new sequence 
of previously defined contour points at all stations of a 
component or segment. Consider a component defined by a 
first $SHAPE set or auxiliary geometric input cards which 
has eleven contour points per station. Assume that it is 
desired to delete point 8 from the contour, add a new contour 
point between 2 and 3, and duplicate contour point 5 in 
order to input a step change in the surface weight at the 
position of contour point 5. The second $SHAPE set for the 
component would contain: 

NPTS = 12, KEEP = 1,2,2,3,4,5,5,6,7,9,10,11, 

The program would construct new section data at all stations 
such that each new station contour would contain the data 
from the original contours taken from the contour locations 
specified by KEEP. The new station 3 may be replaced through 
use of another $SHAPE input set. 

New contour data and subsequent corrections must correspond 
to the revised point numbering. The operation of KEEP is 
bypassed if there is no KEEP input (KEEP(l) = 0) and will 
be ignored if it is input with the first $SHAPE set. Caution : 
The program does not rearrange RHO, H or SURFWT. It is left 
to the user to ensure that these inputs correspond to the 
final numbering system including changes to NSTA. 

KPLT : Integer which controls the mass calculation and con- 

figuration plotting. 

If KPLT = 0, the program computes the mass properties 
and plots the configuration. 

If KPLT =+l, the program plots the vehicle but skips the 
mass properties calculating. 

If KPLT =-l, the program computes the mass properties, 
but skips plotting. 

KPLT does not have to be given the same value for all segments 
but can be varied to skip plotting or mass calculations of 
some segments, etc. The program automatically accounts for 
the number of segments that are to be plotted. The plotted 
segments do not have to be continuous. 

KRHOH : Integer which controls the auxiliary reading of 

elemental density or thickness. If KRHOH < 1, this option 
is skipped. Otherwise, RHO and H input for the current seg- 
ment are expected on formatted cards. The card format is 
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10F7.0 (10 numbers per card; each number can be 7 columns 
wide). H refers to surface elements between the grid points. 
Therefore, the cards contain in order: 

1. NPTS-1 values of RHO for the quadrilaterals between 
stations 1 and 2. 

2. NPTS-1 values of H for the same quadrilaterals. 

3. NSTA-1 repetitious of 1 and 2 for the quadrilaterals 
between subsequent stations. 

Values for H cannot occupy the same cards as RHO and vice 
versa. (See Auxiliary Density/Thickness Input. ) 

KSWT: Integer which controls the auxiliary reading of sur- 

face element weights. If KWST < 1, this option is skipped. 
Otherwise, surface weight distribution (SURFWT) for the 
current segment is expected on formatted cards. The format 
is 10F7.0 (10 numbers per card; each number can be up to 1 
column wide) . The formatted cards contain in order: 

1. NPTS values of SURFWT for the first station contour. 

2. NSTA repetitious of 1 for subsequent station. 

Values of two stations cannot occupy the same card. 


KWNG : An integer which controls the specification of wing 

geometry. It is used in conjunction with the KXYZ option 
described below: 


If KWNG = 0 , the geometry in input in body station 

coordinates in accordance with the KXYZ 
options . 


If KWNG = 1, the upper and lower surface geometry are 
specified in pairs from the leading edge. 
The program converts the input data 
for the pair to a single station. 



INPUT DIRECTION 


INPUT DIRECTION 


If KWNG = 2, The upper surface and the mean line are 
specified with respect to the leading 
edge in pairs. The program computes the 
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the lower surface geometry and replaces 
the mean line input then converts to a 
single station. 

If KWNG = 3, the upper and lower surface are specified 
in a single sequence of contour points. 



Figure 9 summarizes the wing input options available by 
cornering input values of KWNG and KXYZ . 

KXYZ : Integer which controls the auxiliary reading of sur- 

face geometry data for both wing and body stations. The 
following discussion pertains to input of body station contour 
A discussion of wing station contours is given under the he 
heading KNWG. 

If KXYZ < 1, All geometry data is expected in the 

$SHAPE input. There are three auxiliary 
data read options invoked by KXYZ =1, 2 
or 3 respectively. 

If KXYZ = 1, contour inputs for the current segment 

are expected on formatted cards. The for- 
mat is 10F7.0 (10 numbers per card; up to 
7 columns per number is allowed) . The 
formatted cards contain in order: 

1. NSTA values of X for the complete component. 

2. NPTS values of Y for the first station contour. 

3. NPTS values of Z for the first station contour. 

4. NSTA repetitious of 2 and 3 for subsequent stations. 

Values for different names cannot occupy the same card. 

If KXYZ = 2, the program generates contour points for 

a circular cross section from input values 
of radius and vertical offset. 

The contour points are spread in even angular increments 
between the lower center plane and the upper center plane of 
the body. The number of generated points is specified by the 
user (NPTS) . Let I = point number on contour and J = station 
number. The radius of each circular section is input to the 
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KWING KXYZ 


INPUT DESCRIPTIONS 


1 


2 


3 


Pi 

D 

ft 

53 

H 

w 

ft 

<1 

ffl 

CO 

i-w- 


Input upper and lower surface (one station 
for each) from loading edge in $SHAPE 
ZWNG = F (XWNG) . XWNG must be even increments. 


Input upper surface and mean line (one station 
for each) . ZWNG = F (XWNG) for all YWNG. 

XWNG must be even increments. 


Input upper and lower surface as a single 
station. ZWING = F (XWNG) for several 
values of YWING. 


1 


2 


1 


2 


Input wing ordinates in percent chord. 

Input leading edge coordinates and chord length. 
Input camber line in percent chord. 

Input upper surface thicknesses in percent chord. 
Input lower surface in percent chord. 

Input wing ordinates in percent chord. 

Input leading edge coordinates and chord length. 
Input camber line in percent chord. 

Input half thicknesses in percent chord. 

Input wing ordinates in percent chord. 

Input leading edge coordinates and chord length. 
Input upper surface in percent chord. 

Input lower surface in percent chord. 


2 


Input wing ordinates in percent chord. 

Input leading edge coordinates and chord length. 
Input half thicknesses in percent chord. 


FIGURE 9 INPUT REQUIREMENTS FOR WING GEOMETRY. 
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first row of Y and the centerline offset from the X axis 
is input to the first row of Z. Note that this is row , not 
column . The input values are replaced at execution time by 
the appropriate contour point data. 

If KXYZ _> 3, contour points for the current segment are 
expected on formatted cards prepared by a 
digitizer which provides coordinates of 
points (DCP ) directly from drawings. 

The digitizer senses discrete positions of a "reader” placed 
on a drawing. Thus, the digitizer cards contain coordinates 
of points around the station contours in units corresponding 
to the digitizer scale. Additional data is needed to trans- 
form the card input into true units and correct for the 
angular position of the Y,Z axes relative to the digitizer 
axes. The cards must contain in order: 

1. True distance between two points and the DCP of the 
points on the y and z axes. Format (F8. 3 , 4F4 . 0) < 

2. True X coordinate of the first station contour, DCP 
of the Y = Z = 0 point on the drawing, DCP of an 
arbitrary point on the +Y axis, DCP of selected 
points along the contour. Format (F8 . 3 , 16F4 . 0) . 

3. If NPTS > 6, the DCP of the remaining points follows 
on successive cards. Format (18F4.0). 

4. NSTA repetitions of 2 and 3 for subsequent stations. 

Values for different stations cannot occupy the same card. 
Values can start in column 1. Note that each DCP is two 4- 
digit values. (See Digitizer Data Specifications) 

KYZRO : Integer controlling the closure of a section contour. 

Normally, the contour is expected to meet its symmetrical 
counterpart on the Z axis but the input data may not have 
Y = 0.0 exactly at these two points. If KYZRO = 0, the pro- 
gram sets the first and last value of Y at each station = 0.0. 
If the cross section is a closed contour that does not cross 
the Z axis, KYZRO = 1, will set the last contour point equal 
to the first contour point. KYZRO = 2 leaves the input 
unchanged. The KYZRO control is evoked after the input is 
scaled. (See figure 5) 

LABEL (word) : An alphameric label input for each segment 

(and each black box). It can be up to 60 characters long, 
input with a hollerith character count, nH. A long hollerith 
string cannot be continued on a second card, but can be input 
in multiple of 10 character e.g. : 
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LABEL (1) 


20H THIS IS A SAMPLE BLA, 


LABEL (3) = 12H CK BOX LABEL, 

Note that blank spaces are included in the hollerith count. 

NOSE ; An integer which identifies the first station as a 
singular station (i.e. all contour points have the same value) . 

NPTS: The number of points defined around each station con- 
tour in the current segment: '2 £ NPTS £ 45. 

NSEG : The number of components or segments in a complete 

case. It must be supplied in the first $SHAPE set of a case. 

If NSEG = 0, the program skips to processing black boxes. 

NSTA: The number of station planes in the current segment: 

2 £ NSTA <25. 

RFLNTH : The reference length of the vehicle. It is used to 

compute the percent location of the center of gravity c.g. 
and center of volume c.v. The origin for the reference 
length and percent calculation is taken to be the X-coordinate 
of the first station in the first segment. If the preset 
value —100.0 — is used, the percentages become axial distance 
from this origin. 

RHQ (40 , 25) : Density of the quadrilateral surface elements. 

RHO ( I . J . ) is the density of an element defined by four points 
starting at the geometric point (I.J.) and extending one 
point beyond in each direction. 

SCLXs Scale factor applied to all the X-coordinates for the 
current component. It can be used for units conversion as 
well as for changing the vehicle size. SCLX operates on the 
input data after the last $SHAPE set for a segment has been 
read. (See figure 7B) 

SCLY and SCLZ : Scale factors applied to all Y and Z coordinates 

respectively. 

SURFWT (40,25) : Weight per unit area distribution of the 

surface of the configuration. SURFWT is input at every grid 
point. Its contribution to each surface is obtained by linear 
interpolation to determine the mean surface unit weight at 
the centroid of the surface element. 

X (25) : The axial coordinates of the station planes within 

a segment. 


40 



Y (40,25) : The horizontal coordinates of grid points around 

each body station contour. Each column in this array con- 
tains the Y coordinates for all points at one station. Up 
to 40 contour points per station may be used. Up to 25 
stations per component may be specified. 

Z (40,25) ; Similar to Y (above) . Z is the vertical coordinate 
of grid points and can be negative. 

XOF : Offset from the system references axes to the origin 

of the X- values that are input. XOF is added to every value 
of X after scaling (see SCLX) . The result replaces the X 
array. 

YOF and ZOF : Similar to XOF (above) . They are added to 

the Y and Z inputs, respectively. 

XWNG (40,25) : Chordwise coordinates of the grid points around 

each wing-station contour. Each column in this array con- 
tains the X-coordinates for all points at one station. Up 
to 40 contour points per station may be used. Up to 25 
stations per component may be specified. A station may con- 
sist of a complete wing section or only one-half (upper 
surface) and a camber line may be included (see KWNG input) . 

YWNG (40) s Spanwise positions of the wing contours. 

ZWNG (40 , 25) : Vertical coordinates of the grid points around 

each wing- station contour corresponding to the XWNG input 
above. 

Auxiliary Geometry Input . - The program provides for 
auxiliary geometry input following the first $SHAPF, input 
set for each component. This geometry input is basically 
compatible with the geometry input scheme of reference 16 
(Harris input) . Only two of the Harris input options are 
available in VAMP, the fuselage input and wing input. All 
data is placed in the input cards in 10F7.0 format (10 values 
per card and 7 column per input value) . Each new subset of 
data (coordinate points) must begin with a new card. 

Fuselage Data Cards : If the fuselage is circular and cam- 

bered, no auxiliary data cards need be supplied. The radius 
and camber data are supplied in the $SHAPE input set. 

If the fuselage is of arbitrary shape, the y coordinates 
for a half section are given (NPTS values) as Y(I.J.) where 
J is the station number. Following these are the corresponding 
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Z ordinates (NPTS values) for the half section as Z (I . J. ) 
where J is the station number. Each station will have a 
set of Y and Z cards and the convention of ordering the 
ordinates from bottom to top is observed. 


For each fuselage segment a new $SHAPE data set is required. 

If the auxiliary Harris input is specified, a set of cards 
sls described above must be provided. The segment descriptions 
should be given in order of increasing X. 


Wing Data Cards : The first wing data card (or cards) con- 

tain the percent- chord locations at which the ordinates of 
all the wing airfoils are to be specified. There will be 
exactly NPTS percent-chord locations given. (See KWNG Input) 


The next wing data cards (there will be NSTA cards) each con- 
tain four numbers which give the origin and chord length of 
each of the wing airfoils that is to be specified. The cards 
representing the most inboard airfoil are given first, followed 
by the cards for successive airfoils. The information is 
arranged on each card as follows: 

Columns Description 

1-7 x-ordinate of the airfoil leading edge 

8-14 y-ordinate of the airfoil leading edge 


15-21 


z-ordinate of the airfoil leading edge 


22-28 the airfoil streamwise chord length 


If a cambered wing has been specified, (KXYZ = 1) , the next 
set of wing data cards is the mean camber line cards. The 
first card contains up to 10 AZ values, referenced to the 
Z-ordinate of the airfoil leading edge", at each of the specified 
percent of chord locations for the first airfoil. If more 
than 10 values are to be specified for each airfoil, there 
will be NPTS values, the remaining values are continued on 
successive cards. The remaining airfoils are described in 
the same manner, data for each airfoil starting on a new card, 
and the cards arranged in the order which begins with the 
most inboard airfoil and proceeds to the outboard. 


Next are the wing airfoil ordinate cards. The first card con- 
tains up to 10 half-thickness (KWNG = 2) ordinates of the 
first airfoil expressed as a percent of the chord length. 

If more than 10 ordinates are to be specified for each airfoil 
(there will be NPTS values) the remaining ordinates are 


42 



continued on successive cards. The remaining airfoils are 
each described in the same manner, and the cards are arranged 
in the order which begins with the most inboard airfoil and 
proceeds to the outboard. 

If KWNG = 1, the above set of data cards contain the upper 
surface data in precent chord location followed by another 
subset of data containing the lower surface data in percent 
chord locations. (See figure 9) 

Digitizer Data Specification : Fuselage data may be specified 
to the VAMP program through use of data digitizer output. 

The digitizer data option is invoked by specifying KXYZ = 3 
in the $SHAPE input set for the current section. A data 
digitizer is a device which converts graphical information 
into digital data. It measures the relative distance among 
points on a configuration drawing, formats the data and 
records it on an output recorder. When coupled with a key- 
punch machine, the digitizer provides coordinate data on 
punched cards in F4.0 format. The data recorded is suitable 
for computer processing. The VAMP program accepts digitized 
data (KXYZ = 3) and converts it for use in computing mass 
properties of surface elements. 

The data digitizer is usually supplied with a coordinate 
digitizing table upon which the user places the drawing of 
the configuration under study. The following procedure is 
employed in generating digitized geometry input data for 
the VAMP program. 

1. Turn on the power switch. If an absolute zero is 
not provided on the data digitizer, move the curser 
to a convenient point in the lower left portion of 
the digitizer table. Reset the x and y to zero. 

2. Turn on the keypunch machine, set autofeed to on 
and feed two cards down. Depress foot petal to 
test digitizer. Remove the test card from the 
keypunch machine. 

3. Generate a calabration card. This is done by hand 
punching a number representing a convenient measure 
in F8.0 format on the first card as illustrated in 
figure 10. Then digitize two points representing 
that length (measure) on the drawing. 

4. On a new card, hand punch the x-value of the first 
station in F8.3 format. Now proceed to digitize 
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the contour points in a counter clockwise manner 
as illustrated in figure 8. If more than eight 
points are needed, the data will be punched on 
additional cards. 

5. Repeat number 4 above for each station. 

Figure 10 provides an illustration of a data deck digitized 
from a fuselage station. Digitizer data specification can 
not be used for wing stations. 

Auxiliary Density /Thickness Input . - The elemental density 
and thickness are usually provided to the program through 
the $SHAPE input list under the variable names RHO and H 
respectively, but can be read into the program on formatted 
(10F7.0) input cards following the first $SHAPE input as 
illustrated in figure 11. The auxiliary read option is 
invoked by the parameter: 

KRHOH = 1, 

set in the first $SHAPE input list. When the option is 
invoked, input cards must be provided with the following 
information in 10F7.0 format (10 values per card, 7 columns 
per value) : 

1. Values of RHO for the surface elements between 
station 1 and 2. 

2. Values of H for the surface elements between station 
1 and station 2. 

3. Repetition of 1 and 2 for all subsequent stations 
for the current component. 

Each new set of RHO or H must start a new card, (figure 11 ) 

Further, if the auxiliary geometry option described above 
is invoked simultaneously with the auxiliary density/thick- 
nesses input option, the auxiliary density /thickness data 
set follows the auxiliary geometry. 

Auxiliary Surface Weight Input . - The elemental surface weight 
distribution is usually provided to the program through the 
$ SHAPE input list under the variable name SURFWT but it can 
be read into the program on formatted (10F7.0) input cards 
following the first $SHAPE input set similar to the auxiliary 
density thickness option (figure 11). The option is invoked’ 
by the parameter: 
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FIGURE 11 AUXILIARY INPUT SETUP FOR DENSITY/THICKNESS, 


KSWT = 1, 


set in the first $SHAPE input list. When the option is in 
effect, the following information must be provided ten 
values per card and seven columns per value (10F7.0 format): 

1. Values of SURFWT for each point on the geometric 
contour for the first station. 

2. Repetition of 1. for all stations of the current 
component. 

SURFWT values for each new station must start on a new card. 

Any auxiliary geometry and/or auxiliary density/thickness 
data preceed the auxiliary surface weight data in the deck 
setup. The sequence for the three auxiliary data sets is:- 

1. Geometry (KXYZ > 0). 

2. Density/ thickness (KRHOH > 0) . 

3. Surface weight distribution (KSWT > 0) 

Any or all of the above data sets may be omitted assuming the 
appropriate option parameters are set to zero. 

Black Box Data 

The configuration shape data described above provides input 
to the program for analyzing shell structures such as fuse- 
lage, wings, tank, etc. However, the analysis of masses not 
conforming to the surface model or having an unknown shape 
may also be specified by inputting the center of gravity 
(c.g.) and mass properties. 

The mass sources called black boxes may be described by the 
c.g. and weight only (point mass). However, the inertias 
may be taken into account by specifying the moment of 
inertia of the elemented mass about the local axis through 
its own c.g. The orientation of the local axis may be 
specified relation to the system reference or some inter- 
mediate reference system. 

Input Specification for the $BLKBOX Input Set . - The follow- 
ing d i s c uifiT on~l? r o vi d e s a description of the inputs required 
to evaluate mass sources by the black box method. The input 
format is standard FORTRAN namelist as described above (see 
Configuration Shape Data) . 
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The program searches for a $BLKBOX input set after reading 
and analyzing all configuration shape data. The format of 


the data is: 

$BLKBOX 

$BLKBOX -----$ 


The dashes indicate (name=value) sets which will be described 

below: 

BXPLN(3,3): Black box orientation array. This array 

contains the coordinates of three points. 


In matrix notation. 


BXPLN 



X. 




The first column of BXPLN contains the X,Y,Z of the local 
c.g. Second column contains the coordinates of a point on 
the local +x axis. Third column contains the coordinates 
of the local +y axis. Thus, BXPLN defines the position and 
orientation of the black box. BXPLN is used to rotate and 
translate the inertias to the reference axes. If NBX (see 
below) is negative, only the rotation is performed and the 
resultant inertias are saved (see PMIP) and combined with 
the next black box input. Once NBX is set positive, complete 
transformation of the accumulated data is performed. 

CNVMS : Factor to convert weight units into mass units. 

This input parameter affects only the moment and product 
of inertia calculations. The total weight remains in 
input units. The weight elements entered into an inertia 
calculation are divided by CNVMS for conversion to mass. 
The initialization of CNVMS is preset to 32 . 1741b/slug . 


IXX: Mass moment of inertia of a black box about a 

local x axis through its c.g. It must be in consistant 
inertia units. It is not affected by the value of CNVMS. 

IYY and IZZ : Mass moments of inertia for the local y 

and z axes respectively. 
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LABEL ( 6 ) ; An alphameric label input for each segment 
(and each black box). It can be up to 60 characters 
long, input with a hollerith character count, nH. A 
long hollerith string cannot be continued on a second 
card, but can be input in multiple of 10 character c.g. 

LABEL (1) = 2 OH THIS IS A SAMPLE BLC, 

LABEL (3) 12H CK BOX LABEL, 

Note that blank spaces are included in the hollerith 
count. 

NBX: Control integer for black box input. If NBX = 0, 

program proceeds to plotting data or to the next case. 
NBX = 1 is used to specify non- symmetrical transforma- 
tion of the inertias in PMIP to system reference axes 
and add to system total. NBX = 2 is used to specify a 
symmetrical transformation. 

I XY = I YZ = 0,0 

The weight and remaining inertias are doubled, trans- 
lated to system reference axes, and added to system 
total. If NBX is negative, the program rotates the 
inertias but does not translate to the reference axes. 

As long as NBX remains negative, the accumulation of 
inertia properties in PMIP continues. When NBX becomes 
positive, the current contents of PMIP are transformed 
to the system reference axes according to whether NBX 
is 1 or 2. 

PMIP (3,3) ; This array contains the black box' s products 
and moments of inertia about an axis through the local 
c.g. 

In matrix notation, 



T T 

XY YY 
I XZ I YZ I ZZ 


The array is actually symmetrical about the major dia- 
gonal but only the lower terms are input. The program 
fills the upper half. The diagonal terms can be input 
with the names IXX, IYY and IZZ. If the NBX = negative 
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feature is being used, this array is input only with the 
first $BLKBOX set in the group of black boxes being 
accumulated. Input must be in consistent inertia units. 

WTM: Weight of a black box. If NBX = 2, this is the 

weight of the black box contribution to the black box 
set. The mass contribution of WTM to the inertia trans- 
formation is determined by conversion using CNVMS. 

XO : This is the X coordinate of a black box c.g. This 

input will place a value in the X position in BXPLN. 
When NBX is positive, XO must be the black box c.g. 
measured parallel to the absolute system reference axes. 

YO and ZO : Similar to XO (above) . These are Y and 

Z coordinates. Their values are placed in the rirst 
column of BXPLN. If NBX = 1, YO can be negative. 

XOF : X-offset from the system reference to the origin 

of the XO-values that are input, XOF is added to every 
value of XO. 

YOF and ZOF : Similar to XOF (above) . They are added 

to Y and Z inputs respectively. 

Plotting Data 

The program provides optional plotting for all configuration 
geometry on a component by component basis. The integer KPLT 
(set in $ SHAPE) controls the storage of plot data for each 
component. Once the VAMP analysis is complete, any number 
of plots may be generated from the stored geometry for later 
plotting on the CALCOMP plotter (see figure 4) . 

Four plot options are available,, each requiring a single card 
input. The four plot type options and required input for 
each are described below. The program continues to read 
cards and produce plot information until after completing a 
request that included a one (1) in column 72. 

The title card is printed on every plot. The contents of 
the plot request card are printed on all plots except those 
produced by the VU3 option. 

Orthographic Projections . - An orthographic projection 
request card contains the following data punched in the 
indicated columns: 
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Column 

Fortran Name 

Description 

1 

HORZ 

"X," "Y" or "Z" for horizontal axis. 

3 

VERT 

"X," "Y" or "Z" for vertical axis. 

5-7 

TESTl 

The word "OUT" for the deletion of 
hidden lines, otherwise, leave blank. 

8-12 

PHI 

Roll angle, degrees. 

13-17 

THETA 

Pitch angle, degrees. 

18-22 

PSI 

Yaw angle, degrees. 

48-52 

PLOTSZ 

PLOTSZ determines the size of the 
plot. A scale factor is computed 
using PLOTSZ and the maximum dimen- 
sion of the configuration. 

53-55 

TYPE 

The word "ORT." 

72 

KODE 

1 if this is the last plot. 


The orthographic projections are created by rotating each 
point on the body surface to the desired viewing angle and 


then transforming the points into a coordinate system in the 
plane of the paper. The body coordinate system is coincident 
with the fixed system in the plane of the paper when all of 
the rotation angles are zero; for example, the configuration 
X axis and Y axis would coincide with the paper for plots in 
the X , Y paper plane. This would produce a top or plan view 
with the nose to the left. The Y,Z view is from the rear and 
the X, Y view (elevation) is from the Y = negative side, nose 
to left. The rotations of the body and its coordinate system 
to give a desired viewing angle are specified by angles of 
roll, pitch and yaw (0,$,'?) which are positive by the right 
hand rule. HORZ and VERT specify the starting position from 
which the view is rotated, according to the input angles. 

An attempt is made to center the given configuration within 
the specified field. If the desired plot size is greater 
than 28 inches, centering is attempted within 28 inches so 
care must be taken in choosing the view, i.e., PLOTSZ may 
exceed 28 as long as the resultant vertical dimensions of 
the plot do not exceed 28 inches since the paper is 29.5 
inches wide (but up to 120 feet long) . 
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Minimum values are adjusted so that body axis lines with no 
rotation angles coincide with grid lines on the plotter paper. 
Therefore, the plotter pen should be positioned one inch from 
the side of the plotting space and on the intersection of 
heavy grid lines at the start of plotting if plotting is to 
be done on grid paper. The initial pen position is not input 
to the program but specified on the CALCOMP plotting instruc- 
tion card submitted to the computing center after the program 
VAMP run is completed. 

Plan, Elevation and Rear Views (Stacked) . - A three-view plot 
request card contains the following data punched in the 


indicated 

columns : 


Columns 

Fortran Name 

Description 

5-7 

TEST1 

The word "OUT " for the deletion of 
hidden lines, otherwise, leave blank. 

8-12 

PHI 

Y-position on the paper of plan view, 
inches. 

13-17 

THETA 

Y-position on the paper of the side 
view, inches. 

18-22 

PSI 

Y-position on the paper of the rear 
view, inches. 

48-52 

PLOTSZ 

PLOTSZ determines the size of the 
plot. See Orthographic projections. 

53-55 

TYPE 

The word "VU3 . " 

72 

KODE 

If 1, this is the last plot. 


This option produces the X, Y, X, Z and Y,Z orthographic projec- 
tions with no rotation angles. The views are stacked across 
the paper. No attempt is made to center the views or keep 
them from overlapping each other. The final origin point of 
the X,Y,Z arrays is used as the origin of the two axes that 
are being plotted for each view. The " Y-position" defines 
the location of the horizontal axis being plotted for the 
particular view. 

Perspective View . - A perspective view request card contains 
the following data punched in the indicated columns: 
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Columns 

Fortran Name 

Description 

8-12 

PHI 

x of the viewpoint (location of viewer) 
in the data coordinate system. 

13-17 

THETA 

y of the viewpoint in the data coordinate 
system. 

18-22 

PSI 

z of the viewpoint in the data coordinate 
system. 

23-27 

XF 

x of the focalpoint (determines direction 
and focus) in the data coordinate system. 

28-32 

YF 

y of the focalpoint in the data coordinate 
system. 

33-37 

ZF 

z of the focal point in the data coordinate 
system. 

38-42 

DIST 

The distance from the eye to the view- 
ing plane. 

43-47 

FMAG 

The viewing plane magnification factor. 
It controls the size of the projected 
image. 

48-52 

PLOTS Z 

The diameter in inches of the viewing 
plane. DIST and PLOTSZ together 
determine a cone which is the field 
of vision. 

53-55 

TYPE 

The word "PER." 

72 

KODE 

If 1, this is the last plot. 


Perspective views are obtained by projecting straight lines 
from points on the object to the viewer; the viewer being a 
single point in space. The points where these lines intercept 
the view plane form the image. If the view plane is between 
the viewer and object, the image is smaller than the object 
because the projected lines are converging. Thus, the basic 
scale of the plot is controlled by the ratio of DIST to the 
viewer's distance from the object. In this program, the view 
plane cannot be placed behind the object. It can be placed 
within the object, in which case, a cutaway view is obtained 
of what remains behind the plane. The focal point is the 
point towards which the viewer is looking; also called center 
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of view. The view plane is normal to the line from the viewer 
to the focus. PLOTSZ is the diameter of a circular "hole" 
located at DIST, through which the object is viewed. Per- 
spective line intercepts outside this circle are not plotted. 
Thus, PLOTSZ is the maximum size of the actual plot, but it 
can be smaller. FMAG magnifies the view of the object (FMAG 
can be less than 1.0). It affects the dimensions normal to 
the viewer-focus line. The magnified view is still seen 
through the PLOTSZ hole. If DIST and PLOTSZ are such that 
the entire object cannot be seen, even less will be seen as 
FMAG increases. Each plot will be inside a PLOTSZ circle. 

Data which are within the cone of the viewing plane, but 
not in the immediate range of the focal point, may be dis- 
torted. 

Stero Frames Suitable for Viewing in a Steroscope . - The input 
is identical to that for the perspective views except that 
the word "STE" is used in columns 53-55. DIST, FMAG and 
PLOTSZ must be selected to produce a plot that will fit the 
intended stereo viewer. The stereo plots are two perspective 
views with the viewpoint shifted. These two plots are 
centered 9.8 inches apart. 

Use of VAMP within ODIN 

The Optimal Design Integration (ODIN) system is a library of 
independent computer programs representing the analytical 
capabilities in a wide variety of technological disciplines. 

The VAMP computer program is but a single member of the ODIN 
library. The sequence of execution of the individual pro- 
grams is controlled by the executive program, DIALOG (reference 
15) which also maintains a name-oriented data base of design 
information. Each piece of information is stored by name. 

The data base forms a communication link among the programs 
in the library. When used within the ODIN system, VAMP 
receives data from the data base before execution and pro- 
vides information for updating the data base. The actual 
transfer of' information to and from the data base is per- 
formed by DIALOG through pre- and post- processing of the 
data so VAMP is "unaware" that it is part of an analysis 
involving many programs. The general procedure for using 
ODIN is described in reference 15, but there are no special 
input requirements for using VAMP within the ODIN system. 

A single control directive, 

’EXECUTE VAMP’ 

is required for the execution of the program. The delimiter 
C) is a 4-8 punch. The data which follows this directive 
is the normal input data described above. However, any data 
values may come from the data base by specifying the data 
base name on the input card (in lieu of the actual value) . 


54 



$SHAPE 

I 

I 


I 


SCLX = ' SCALE , ' 
SCLY = ’ SCALE , ! 
SCALX = ' SCALE , ’ 


$END 

In the above illustration the name SCALE is a data base name 
which may represent a scale factor for the entire configura- 
tion. The executive program,, DIALOG replaces the name and 
the associated delimiter 'SCALE' with the data base value 
for SCALE. Upon execution of VAMP, the input component is 
photometrically scaled by the current scale factor in the 
data base. A procedure is also available for transferring 
data base arrays (see reference 15) . 

The illustration above applies to namelist input, but the 
procedure for extracting data base information is equally 
applicable to formatted input. The field width is specified 
to the position of the delimiters (') as illustrated below: 

/ 2579 3.426 'DEN' 7.29 


DEN is assumed to be the name of a data base variable. The 
value of DEN will be placed (by DIALOG) on the input card in 
the most significant (E or F) format left justified in the 
specified field. If the data base variables were an integer 
the value would be right justified in the field (I format) . 
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PROGRAM OUTPUT 


The output from the VAMP program is divided into five 
groups . 

Input deck listing 
Configuration data output 
Black box data output 
ODIN output 

Configuration plot data output 

The first three groups are printed in the listed order on 
the normal output elements. The ODIN output is a subset 
of the configuration and black box data saved or a special 
output file for potential storage in the ODIN data base. 
Configuration plot data are stored as a plot vector file 
(a physical tape) for later plotting on a CALCOMP device. 

A sample case illustrating all data output is presented in 
Appendix A. 


Configuration Data Output 

The program computes the volume, areas and mass properties 
of shell structures on a section-by-section basis for each 
component and on a segment- by- segment (component) basis for 
the entire configuration. It accumulates the results of 
previous (sections, or segments) with the current section 
(or segment) as the calculation proceeds. One page of out- 
put presents the mass properties of the current section and 
the cumulative results of all previous sections . 

Intermediate results for each component are summarized by 
a page of output presenting the mass properties of the 
current segment and the cumulative results for all previous 
segments. The last such summary provides mass properties 
for the entire configuration. The mass properties output 
for each segment is preceeded by a printout of program con- 
trols and corner-point geometry. The first page contains 
the current values of all the program option controls, the 
X coordinates of the segment, the cross section areas of all 
station planes and the net half widths of all station con- 
tours. A new page is started with the corner-point geometry 
of the segment. If mass property calculations were requested 
(KPLT<0) , H, RHO and SURFWT are then printed for the segment. 
The geometry and unit mass data is printed in station order 
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(J) down the page and around the contour (I) order across 
the page. The format provides for 11 points to be printed 
across the page. Since the program can accept up to 45 
points per contour, the additional points follow in separate 
groups. The (I,J) indices for each value's position in its 
respective array are clearly indicated so that the proper 
subscripts can be readily determined if changes in the in- 
put are needed. 


Black Box Data Output 

The black box data output is similar in structure to the 
configuration data output. The elemental properties of the 
black boxes are rotated and/or transformed to the system 
references one at a time until all black box input is satis- 
field. The program prints the individual properties and 
the accumulated properties as the calculation proceeds. 

The accumulated properties include all configuration data 
as well as the properties of previous black boxes. The 
accumulated properties for the last black box represents 
the mass properties for the entire vehicle. 

The black box output contains values printed in columns. 

The columns are identified at the top of each page. The 
columns include X,Y,Z coordinates, weight and the six 
moments and products of inertia. Each line starts with an 
alphabetic description. The second output line starts with 
the word INPUT and presents part of the input supplied for 
the current black box. The X,Y,Z columns are the coordinates 
of its c.g. (the first column of BXPLN) . The inertias in 
this line are the input values, relative to the local prime 
axes. The third line starts with: 

ROT- TRANS /REF . AX 

This line contains the mass properties of the black box 
about the system reference axes. If no rotation of the 
inertias was required, the letters ROT are deleted. If 
NBX = 2, the current black box has a symmetrical counter- 
part and Y = 0.0 on this line along with the appropriate 
doubling or zeroing of the weight and inertias. After trans- 
forming the black box to the reference axes, program VAMP 
adds it to the cumulative system which includes all segments 
of the vehicle and all proceding black boxes. The resultant 
mass properties about the reference axes, are printed on 
the line that starts with SYSTEM/REF. AXES and continues with 
no entries in the X,Y,Z columns. The fifth (and last) line 
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leads off with SYSTEM/C. G. and contains the c.g. coordinates 
of the accumulated system, total weight and inertias about 
axes through the system c.g. , parallel to the reference axes. 

If NBX is negative, only the first and second lines are 
printed (LABEL and INPUT) but the program expects more input 
to complete the current black box. With each succeeding 
$BLKBQX set, the INPUT line reflects the rotation specified 
by the preceding set. ROT is still deleted if the final 
$BLKBOX set does not specify rotation. The five line for- 
mat is repeated for all black boxes in the case. 

Configuration Plot Data Output 

Plotted output from the VAMP program depends on the occurrence 
of the following four events: 

1. The setting of KPLT >_ 0 for at least one of the 
configuration components. 

2 . The presence of at least one plot request card 
following last black box data input set. 

3 . The transfer of the plot vector file CALTPE to a 
physical tape as illustrated in figure 4. 

4 . The submittal of a CALCOMP plotting instruction 
card. 

Any of the four types of plots described earlier may be 
generated. Figures 1 and 2 illustrate CALCOMP output for 
the program. Appendix A also contains a sequence of CALCOMP 
output examples for the sample case. 

In addition to the plotted data, the program provides some 
printed information also. 

The program prints the minimum, maximum and difference 
between maximum and minimum values of X,Y and Z respectively. 
There are the values used to scale the plots. The program 
also lists the plot request cards. Each card is listed 
twice. First is the card image which shows what was punched 
on the card. The second listing is the computer's inter- 
pretation of the card. This interpretation is designed to 
indicate the correct columns for the data. Thus, it shows 
input rounded off to a whole number, but all decimal places 
are retained internally. It also shows the signs which 
have been known to be punched in the wrong column and lost. 
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If plotting trouble is encountered, the last request card 
listed is the one the program is trying to satisfy. 

Remember that the last plot request card must have a 1 in 
column 72 unless it is also the last card in the input 
deck. If there are no plot request cards and no following 
case data (i.e. end of input file) , the program will make 
a normal stop after printing the maximum/minimum values. 

The actual plots are made on a CALCOMP plotter. Each plot 
provides an image of the title card. With the exception of 
plots produced by the VU3 option, each plot also contains 
the image of the plot request card. Thirty inch plot paper 
should be requested. 


ODIN Output 


The ODIN data base output consists of the volume, area and 
mass properties of each segment on an individual and 
cumulative basis. The information is placed on a special 
file in a namelist like format for use by the DIALOG 
(reference 15) executive program. Any subset of the infor- 
mation in this file may be placed in the data base by 
DIALOG for use by other programs. A special name generating 
capability is employed in VAMP to provide unique descriptions 
(names) for the various segment, component and black box 
properties. 

A set of array information is generated for each segment 
as follows: 


SEGTOTOn 
SEGXYZ On 
SEGSBTOn 

CUMTOTOn 

CUMXYZOn 

CUMSBTOn 


Individual Segment Information 


Cumulative Information 


The character position, n, varies as the segment number 
changes. There are as many of the above sets as there are 
segments in the mass properties analysis. Figures 12 and 
13 associate the elements of the above arrays with the 
actual output from VAMP. The reader unfamiliar with the 
illustrated output is referred to Appendix A for complete 
familiarization. 
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FIGURE 12 ODIN DATA BASE OUTPUT FOR SEGMENT 1 -—FUSELAGE. 



CSGT 001 HARRIS TYPE INPUT OCT 11, 1972 
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FIGURE 13 ODIN DATA BASE OUTPUT FOR SEGMENT 2 WING 




Figure 12 is the normal output for SEGMENT 1 FUSELAGE. 

Generally the upper half of the page is segment data and 
the lower half is cumulative (all previous segment calcula- 
tions) data. Since figure 12 is the first segment, the 
lower half data is identical with the upper half data. The 
generated ODIN data base array names appear on the left in 
bold type and the element numbers of the generated names 
which are associated with the^data elements appear directly 
beneath the corresponding data elements. For example, the 
c.g. location of segment 1 (SEGXYZ01) is: 

X cg = 1.838 = SEGXYZ 01 (1) 

Y =0. = SEGXYZ 01 (2) 

eg 

Z = 1. 966 = SEGXYZ 01 (3) 
eg 

The corresponding cumulative c.g. location CUMXYZ01 is 
illustrated in the lower half of figure 12. Similar 
examples may be observed for surface areas, inertias, etc. 

Figure 13 illustrates the correspondence between VAMP data 
and ODIN data base elements for the second segment. Here 
the upper half is the calculated properties for the wing 
segment. The lower half is the accumulated properties for 
the wing and fuselage (segments 1 and 2) . The ODIN data 
base names are similar to those in figure 12 differing 
only in the last character positions. 

n = 2 for the wing 

n = 1 for the fuselage 

The analyst may quickly determine from the segment number 
printed what ODIN data base name will be generated. Once 
determined the information may be defined in the ODIN data 
base for use by other programs. Up to 99 segment names may 
be generated (i.e. SEGXYZ 09 9) . 

In addition to analyzing configuration segments, VAMP gen- 
erates mass properties (individual and cumulative) for 
elemental point masses (black boxes) . Figure 14 illustrates 
ODIN data base output for one black box. 

The first set of print statements represents the cumulation 
of all properties for all segments. This data is identical 
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FIGURE 14 OWN DATA BASE OUlPUi FOR ONE (1) BLACK 




to the data accumulation for the last segment (see figure 
12) . The ODIN data base output names for this cumulation 
will always be: 

CUMTOTOO 

CUMXYZOO 

with 00 as the last two character positions. Therefore, 
even though the number of segments may change from run to 
run, the vehicle totals (without black box) will always 
bear the above names. 

The second set of data, in figure 14, is for the black box, 
both individually and cumulative. The correspondence of 
the ODIN names and data elements with the black box data is 
indicated below. The individual black box properties have 
root names. 

BBXXYZnO 

BBXSBTnO 

The black box number is identified by n. For example, the 
third black box would be: 

BBXXYZ30 

BBXSBT30 

The cumulative properties of the black boxes have root names 

of: 


CUMTOTnO 

CUMXYZnO 

CUMSBTnO 

The root names are identical to the root names for the seg- 
ment cumulative properties. The differences lies in the 
last two character positions: 

CUMXYZ05 = Cumulative c.g. properties for the fifth 
segment. 

CUMXYZ50 = Cumulative properties of all segments plus 
five black boxes. 
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Up to 99 black box names may be generated. The last 
character position will always be zero (0) such test the 
99th ODIN name would be: 

CUMXYZ990 

The ODIN data base output scheme described above is adequate 
for extracting most any combination of vehicle properties. 
Intermediate or section data is not available directly but 
can be obtained simply by grouping sections into smaller 
segments. The new segments would consist of portions of an 
original segment. For example, the fuselage may be split 
into multiple segments by input. By judicious selection of 
segments, the analyst may obtain the properties desired for 
storing in the ODIN data base. 
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CONCLUDING REMARKS 


Program VAMP is designed to bridge the gap between the 
historical approach to estimating weight and the detailed 
structural design approach. It computes the mass proper- 
ties, center of gravity location, enclosed volume, wetted 
area and planform area of any closed structure that has a 
plane of symmetry. The configuration is described to the 
program by ordered sets of coordinates on its surface. 

The surface is approximated by quadrilateral surface 
elements among the points. The mass properties of the con- 
figuration are computed from properties input at each 
point. The computed mass property totals contain the con- 
tribution from the distributed mass in the vehicle surface 
wall. Additional masses may be added by specifying each 
center of gravity location and mass properties. Picture- 
like images of the configuration under study may be obtain- 
ed from the surface geometry. 

The program can be interfaced with other technology modules 
within the ODIN system. VAMP itself is not a sizing or 
structural synthesis program; other modules in the ODIN 
program library carry out these functions. As a contri- 
butor to the overall design analysis, VAMP provides a 
means for assembling shape and mass properties data from 
other sources and for combining the information to pro- 
duce the total mass properties of a vehicle design. 


The significant contribution of the program to the general 
problem of computer aided design is twofold. Eirst, the 
program offers a unified approach to the accounting of 
the vehicle component mass properties. The fact that a com- 
ponent weight is not well known does not exclude it from 
being analyzed in a mass properties evaluation. Furthermore 
detailed mass properties of some vehicle components may be 
determined elsewhere and still be included in the VAMP 
analysis. Second, the distributed mass surface model 
incorporated in VAMP allows a very convenient means of pro- 
ducing vehicle geometric perturbations early in the design 
process with some degree of confidence in the resulting 
affect of the perturbation. 
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APPENDIX A 


SAMPLE PROBLEM 


The configuration being analyzed with the sample problem is 
a wing body configuration representing an advanced trans- 
portation concept. Both the wing and the body are described 
to the program through use of the auxiliary geometric data 
input option (KXYZ > 0) . The wing geometry is described to 
the program using the camber and thickness distribution 
option (KWNG = 2, KWYZ = 2). The surface weight option is 
employed for specifying the unit weight distribution for 
both the fuselage and wing surfaces. Surface weight dis- 
tribution is described to the program through the $ SHAPE 
input set (KSWT = 0) . 

One black box is input representing the weight of 10000 
pounds located 160 feet aft of the nose and ten feet above 
the fuselage reference line. Five plots are requested. 

The first one is a three-view and the other four are ortho- 
graphic projections representing different views and using 
different hidden line options. The data presented is a 
sequential listing of the data output from the computer. 

The ODIN output is appended by copying the special output 
file to the printer at the end of the job. The geometric 
pictures were obtained by submitting the plot vector tape 
file to the CALCOMP plotter. These pictures are presented 
at the end of the section. 
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